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Abstract 



It has been proposed that the large cross sections for prompt ^ip, ip' , and 
Xc production at the Fermilab Tevatron pp collider can be explained by a 
dominant color-octet term in the fragmentation function for a gluon to split 
into quarkonium. We show that this mechanism makes testable predictions for 
double-quarkonium ip^, tp^' , ipXc ip'^ and TpXb production, as well as for 
production, using color-octet matrix elements previously determined from 
charmonium production data. The ipip signal would already be measurable at 
the Tevatron, while the ipXc and Wip signals would be on the edge of present 
detectability. 
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Typeset using REVT^ 



Until recently the standard way to calculate ip production was the color-singlet model [|l| , 
where the is treated as a cc pair in a color-singlet ^5*1 state with vanishing relative velocity. 
Although the color-singlet model has been sucessful in some applications, it predicts rates 
that fall orders of magnitude below the data when applied to prompt production at 
the Fermilab Tevatron via the hard scattering subprocess g + g ip + g. A proposal to 
explain this discrepancy is that the dominant contribution to charmonium production at 
high transverse momentum px comes from hard gluon production followed by fragmentation 
of the gluon to a cc pair in a color-octet state [P,^. The conversion of the cc pair to ip is 
nonperturbative and involves the absorption or emission of two soft gluons. Because there 
are new, nonperturbative parameters associated with the color-octet mechanism p,^^, the 
rate for single quarkonium production cannot be immediately predicted, but rather the 
experimental data can be used to determine these parameters. In the present paper we 
show that the color-octet mechanism now makes testable predictions for double-quarkonium 
ipip, ipip', ipXc, and ipXb production, as well as for Wip production at the Tevatron (CM 
energy y/s = 1.8 TeV). 

Using the factorization approach developed in Q the fragmentation function for a gluon 
to split into a QQ quarkonium state H is 

D,^H{z,fi')=Y.d9^n{z,fi'){0^) , (1) 
n 

where z = {Eh +PH)/{Eg + Pg) is the light-cone fraction, fi = 2mQ is the renormaliza- 
tion scale, while n denotes the angular-momentum quantum numbers "^^^^Lj and the color 
quantum number 1 or 8. The nonperturbative matrix elements (O^) represent the inclusive 
probability for forming the state H from the QQ pair [Q; color-singlet values may be cal- 
culated from potential models but color-octet values are unknown. The relative importance 
of the different matrix elements may however be determined by how they scale with v, the 
typical relative velocity of the heavy quarks in H. The short distance coefficients dg^n can 
be calculated using perturbation theory in a^; they depend on z, /i and the quantum num- 
bers indexed by n. Note that the dependence on the quarkonium state H appears only in 



the factor (C»f ). 

The leading-order perturbative calculation of the color-octet short-distance coefficient 
gives 

dg^s^s, = -z)^ K{g ^ 8'S^)6{1 - z) , (2) 

with numerical values K{g 8^5*1) = 0.01 GeV^^ for charmonium and K{g 8^Si) = 
2 X 10""^ GeV~^ for bottomonium. Here we put a;s(4m^) = 0.26 and as(4m5) = 0.17, 
with rric = 1.48 GeV and = 4.88 GeV; these masses are taken from the BuchmiiUer- 
Tye potential model calculations of Ref. [0. The leading-order color-singlet term in the 
fragmentation function, which does not have a simple analytic form is of order and 
gives much smaller factors K{g l^Si) = 4.3 x 10"^ GeV~^ (3.5 x 10~* GeV"^) for 
charmonium (bottomonium) cases. It is consequently assumed that the color-octet term 
in the fragmentation function dominates over the color-singlet term, even though (Ot) IS 
expected to be smaller than {Of) by a factor of order v^. 

Numerical values of the nonperturbative matrix elements are given in Table I for the 
S'-wave states ip, ip' and T, and in Table II for the P-wave states Xc and Xb- The color- 
octet matrix elements {Oq{^Si)), {Oq {^Si)) and {Og"'-' {^Si)) are determined empirically 
from Tevatron data |P,|^J^ , based on calculations using both color-octet and predicted color- 
singlet contributions in lowest order. The bottomonium matrix elements {Oj^^^\^Si)), 
{Oj^'^^\^Si)) and {Og''''^^^\^Si)) were estimated from the charmonium values by scaling 
arguments . All our color-octet values follow Ref. (but we have corrected a misprint in 
the value of {0^"^^^\^Si))). The color -singlet matrix elements, shown here for comparison, 
are determined from potential model calculations using a Buchmiiller-Tye potential 0; they 
are related to the radial wave function R and its derivative at the origin by 

(Of(35i)) = (2J + l)^|i?(0)r, i^ = ^,^',T, (3) 
(0f(3p,)) = (2J+l)-^|i?'(O)|^ H = Xc,Xb, (4) 

ZTT 



where A^c = 3 is the number of colors. 

Given these results, predictions can immediately be made for other processes involving 
hard gluon production with the gluon fragmenting to quarkonium. Two recent papers have 



considered prompt Tp and T production at the LEP e+e~ collider via Z qqg decays P,p!0 
In this paper we address tests of the octet fragmentation mechanism that can be made at the 
Tevatron pp collider, based on multiple quarkonium production or on quarkonium production 
in association with a VT-boson. 

In our calculations we shall neglect the evolution with /i^ of the gluon fragmentation 
function, which greatly simplifies the calculation. This introduces some error, but including 
evolution would not necessarily be an improvement, since naive Altarelli-Parisi evolution 
does not respect the phase space constraint Dg^^{z, fj,"^) = for z < M^/fi'^ [|ll| . 

Double-^/? production 

The pt spectrum of prompt single production is presumed to be dominated by the 
color-octet term in the g fragmentation function given in Eq. (|l|). Assuming that 
QQi 99 ~^ 99 s-re the dominant contributing hard subprocesses, we can predict double-'?/' 
production, with both gluons fragmenting to ip, through the ratio 

dpT{Vi)dpT{V2) 2 dpTiv) 
in the absence of cuts, where the | is a combinatorial factor. This relation remains true in 
the presence of pt and angle cuts, provided they are applied equally to g and ip- 

We note that the gq gq and gq gq subprocesses are not in fact completely negligible 
and contribute about 20% of single-'?/' production with pri'P) > 4 GeV at Tevatron energies; 
this correction reduces our prediction in Eq. (^) by about 20% and we take it into account 
below. We note also that the differential cross section on the right-hand side of Eq. refers 
to all single prompt-?/) production (including both direct g ^ ip fragmentation and indirect 
g {XcJi'P') "P contributions but excluding -p^s from the decay of 5-mesons); however, 
the second 'p on the left-hand side is produced directly (excluding Xcj, 'P' and B-meson de- 
cays). To obtain the total prompt double- ?/> rate, we must add similar contributions from Xcj 



and ip' production (see Table III) multiplied by the corresponding Xcj — ^ V'T s-nd ip' ~^ i'^ 
branching ratios, that give approximately the same Prii^) distributions and increase the 
total rate by a further factor ^2.0. In practice is usually detected via ip — > fj^^fi^ decay, 
with branching fraction B{iIj = 0.0597(25); for this decay mode (indicated by the 

notation V'M/i) we obtain the differential cross section for total prompt double-'?/' production 
as a function of pri'ip) shown in Figure 1. The curve was generated from qq, gg — > gg subpro- 



cesses using the MRSDO parton distribution functions [|T2| with renormalization scale and 
factorization scale both chosen equal to the transverse momentum of the fragmenting gluon 
Prig) — Pri'ip)', it includes the enhancements from indirect contributions. A pseudorapidity 
cut of |?7| < 0.6 was imposed on the ip^s produced. 

Integrating Eq. (^ and including the factor 2 enhancement from indirect g {XcJy '^') ~^ 
ip contributions, we obtain 

^{PP^ M,,X) ^ ^ 8'S^){OtCS,))B{^ /iV) = r.5 X lo-*^ , (6) 
a{pp ^i^i^gX) 

where the cross sections are defined with a minimum p^ requirement on all ip. This equation 
remains true after pseudorapidity cuts, provided the same cuts are applied equally to ip and 
to the recoil gluon jet g m pp ^ '^'Q^- The observed cross section for single prompt %p 
production at the Tevatron with pri"^) > 4 GeV and \ri\ < 0.6 is (7{pp —>■ ip^^gX) 24 nb 
0, from which we infer that the cross section arising from gg final states is about 19 nb and 
hence a{pp — > ip^^ip^^X) ^ 0.14 pb. Thus Eq. (|D predicts that of order 10-20 douh\e-ip^^ 
production events should be detectable already with pri^P) > 4 GeV and \ri{ip)\ < 0.6, for 
integrated luminosity of order 100 pb"^ now accumulated by each of the Tevatron detectors. 
We note that experimental factors such as detector efficiency could reduce this number 
considerably. On the other hand, ii ip ^ ee decays are also detected (there should be good 
efficiency for this in events that are already tagged by one ip —>■ fifi), the event rate will be 
increased by a factor 3 or 4, depending on whether double-'?/'ee final states are included. 

tp -\- quarkonium production 

All the dependence on the ip state in Eq.(H) appears only in the matrix element {Og (^Si)) . 



It is therefore trivial to generalize this equation to the case where one gluon fragments to 
a if) while the other gluon fragments to some other quarkonium state. We have tabulated 
the results in Table III. Note that the Xcj and XbJ detection modes and branching fractions 
remain unspecified, and that the g —>■ H fragmentation picture is applicable only for prig) > 
mil- The other branching fractions used are: B{ip' f^^fJ'^) = 0.0077(17), B(T{1S) — >■ 
^+^-) = 0.0248(07), and B{T{2S) fx+fi-) = 0.0131(21). The ip + {xc ^ ipl) channels 
account for about 33% of the prompt double-'?/' production rate, as already noted above, so 
they predict about 5 ip^^ip^^'y events for 100 pb~^ luminosity; there is some extra loss of 
efficiency from the need to measure the photon and reconstruct the Xc invariant mass, but 
these signals would appear to be approaching detectabihty. Similarly, the ij+i^ip' ~^ ipTiir, iprj) 
channels account for about 17% of the prompt double-'?/' production rate, with a few events 
at present luminosities, but it may well be impracticable to reconstruct ip' in these modes 
because of backgrounds. The predicted rates for i'fj,f_iip'^^, fi^T^ fj,^ and ip^^Xb production 
are all too small to be observable at present. However, we note that the indirect process 
g — ^ Xfej(l-P) ~^ T(IS') is predicted to dominate by a factor 3 over direct g T(IS') 
fragmentation, when the relevant Xbj{^P) ^(l'S')7 branching fractions are folded in, so 
the net T(IS') rate is 4 times larger than the third row of Table III would suggest. 

+ gauge boson production 

Another test of the dominance of the color-octet term in the gluon fragmentation function 
is high-pr + gauge boson production. Requiring the ip and gauge boson to have large 
relative pt severely restricts the number of possible production mechanisms. 

High-pT ip + 'J production at the Tevatron, via gluon-fusion and color-singlet- 



fragmentation mechanisms, was investigated in Ref . [|T^ ; it was found that gluon-fusion was 
the dominant production mechanism for all px accessible to experiment. We have investi- 
gated the additional contribution from the subprocess qq ^ 'yg, with the gluon fragmenting 
to ip via the color-octet mechanism. We find that this color-octet contribution falls well 
below the gluon-fusion contribution for pj- < 20 GeV, essentially because there are many 
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more gluons than quarks in the proton at the relevant x-values for Tevatron energies. Hence 
high-pr V + 7 production is not sensitive to the color-octet contribution, at values where 
most events will be seen. 

In contrast, high-p^ production oi if) + W offers much cleaner tests of the color-octet 
mechanism. It has no gluon-fusion contribution, and the dominant direct production sub- 
process is qq Wg with the gluon fragmenting to ip via the color-octet mechanism. This 
mechanism contributes the differential cross section 

^^.(pp ^ Wi;,,X) ^ K{g ^ S,){Ot{' S,))B{i, ^ ^^^^^-)^(pp ^ WgX) , (7) 



dpT dpT 
where pt refers on the left to and on the right to g, for sufficiently large px- Rapidity 
cuts must be applied equally to ip and to the recoil gluon jet. Equation (|^ describes only 
the direct g ^ ip contribution; indirect g {Xcj-,'4'') ~^ 4' contributions will enhance the 
rate by a factor 2 as discussed above. Integrating Eq. (|^ for px > 5 GeV and including this 
enhancement factor, we obtain 

a{pp iljf,f,WeuX;pT{ip) > 5 GeV) ^ 

2K{g ^ 8^S,){OtCS,))B{^ ^ f^^f^'Mpp ^ WeugX-prig) > 5 GeV) ^ 10 fb , (8) 

where both —>■ fifi and W ^ ev branching fractions are included. The value for a{pp — *■ 
WeuO', Prig) > 5 GeV) on the right-hand-side of Eq. cannot be taken directly from 
experiment, since VT+jet production has substantial contributions from gq{q) — Wq'{q') 
subprocesses; it was calculated at tree level for y/s = 1.8 TeV (divide the rate by 2(4) if 
Pt > 10(15) GeV instead). This cross section implies one event per 100 pb^^ and is therefore 
on the edge of detectability at the Tevatron at present; however, including the ip —>■ ee and 
W ^ fiu decay modes would increase the cross section by a factor 4, and future increases 
in luminosity could make it measurable. 

Summary 

The proposal that prompt charmonium production at the Tevatron is dominated by a 
color-octet term in the gluon fragmentation function introduces new parameters that 



can however be determined from data @J^,§]. We have shown that this mechanism can now 
be tested through its predictions for double-charmonium and charmonium+VT production. 
It is remarkable that indirect g — > [Xcj,"^') —>■ contributions enhance signals by about 
a factor 2 and analogous g XbJ contributions, inferred by scaling arguments 0, 

enhance T(IS') production by a factor 4. We have found that the predicted rate for ipi/) 
production should already be observable at the Tevatron, with the pt distribution shown in 
Fig.l; the rates for ipip', ipXc, 4''^^ ^ind ipXb production (summarized in Table III) are smaller 
but the + (Xc i'l) signal could already be detectable too. We have also investigated ip 
plus gauge boson production at high px] + 7 production proves to be insensitive to color- 
octet contributions, but ip + W production offers a clean test of the color-octet mechanism, 
with a predicted rate on the edge of present observabihty at the Tevatron. 

Acknowledgments 

RJNP thanks Paul Dauncey for helpful information and the University of Wisconsin for 
hospitality during part of this work. VB thanks the Institute for Theoretical Physics at 
the University of California, Santa Barbara for hospitality during the completion of this 
work. This research was supported in part by the U.S. Department of Energy under Grant 
No. DE-FG02-95ER40896, in part by the University of Wisconsin Research Committee with 
funds granted by the Wisconsin Alumni Research Foundation, and in part by the National 
Science Foundation under Grant No. PHY94-07194. 



8 



REFERENCES 



[1] J. P. Leveille, Proc Topical Workshop on Production of New Particles (V. Barger, 

F. Halzen, eds), Madison 1979; W.Y. Keung, Proc Cornell Z° Workshop 1981; 
E.L. Berger and D. Jones, Phys. Rev. D23, 1521 (1981); R. Baier and R. Ruckl, Phys. 
Lett. 102B, 364 (1981); L. Clavelh et al, Phys. Rev. D29, 57 (1984); for a review, see 

G. A. Schuler, CERN-TH.7170/94. 

[2] CDF Collaboration: F. Abe et al , Phys. Rev. Lett. 69, 3704 (1992); 71, 2537 (1993); 
K. Byrum, Proc. XXVII Int. Conf. on High Energy Physics, Glasgow, Scotland, July 
1994 (P.J. Bussey, I.G. Knowles eds, lOP Bristol) p.989; FERMILAB-CONF-94/136-E 
( |hep-ex/9412013| ). 

[3] E. Braaten and S. Fleming, Phys. Rev. Lett. 74, 3327 (1995). 

[4] G.T. Bodwin, E. Braaten, and G.P. Lepage, Phys. Rev. D51, 1125 (1995). 

[5] P. Cho and A. Leibovich, [hep-ph/95053^. 



[6] 
[7] 



[9 

[lo: 
[11 

[12 
[13 



M. Cacciari, M. Greco, M.L. Mangano and A. PetreUi, Phys.Lett. B356, 553 (1995). 
E. Eichten and C. Quigg, Phys. Rev. D52, 1726 (1995). 
E. Braaten and T.C. Yuan, [hep-ph/9507398 . 



K. Cheung, W.-Y. Keung and T.C. Yuan, |hep-ph/9509309 
P. Cho, |hep-ph/9509355 . 



E. Braaten, M. Doncheski, S. Fleming and M.L. Mangano, Phys. Lett. B333, 548 (1994). 
A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Lett. B306, 145 (1993). 
D.P. Roy and K. Sridhar, Phys. Lett. B341 413 (1995). 



9 



TABLES 

TABLE I. Values of the nonperturbative matrix elements used to calculate ■0, if)' and T produc- 
tion at the Tevatron. The color-octet matrix elements for charmonium were determined empirically 
from CDF data; the color-octet matrix elements for bottomonium were determined by re-scaling 
the corresponding charmonium matrix elements [5]. The color-singlet matrix elements, shown here 
for comparison, were determined from potential model calculations [7]. 







v' 


T(15) 


T(25) 




1.2 GeV^ 


0.75 GeV^ 


9.3 GeV^ 


4.5 GeV^ 


(Ose^i)) 


0.012 GeV^ 


0.0073 GeV^ 


0.01 GeV^ 


0.006 GeV^ 



TABLE II. Values of the nonperturbative matrix elements used to calculate Xc and Xb pro- 
duction at the Tevatron; different J-values differ simply through a factor (2J-I-1). The color-octet 
matrix element for Xc was determined empirically from CDF data; the color-octet matrix element 
for Xh{^P) was found by re-scaling the Xc matrix element [5]. The color-singlet matrix elements, 
shown for comparison, were determined from potential model calculations [7]. 



Matrix Element 


XcJ XbjO-P) 


(OiePo))/(2J + l) 


0.11 GeV^ 2.0 GeV^ 


{Os{^Si))/{2J + l) 


0.005 GeV^ 0.014 GeV^ 


TABLE III. The ratio of the 


cross section for prompt ip + H production to the cross section for 


prompt V'+5 production with H ■ 


= ijj', XcJ, T(IS'), T{2S),XbJ, based on direct g ^ H fragmentation 


contributions only. 




H 


a{pp ilj^^HX)/a{pp tpij.,j,gX) 




K{g ^ 83^1) B{^' ^ ^l+^l~){Ot'eS,)) = 6 x 10"^ 


XcJ 


K{g 8^ Si) = (2J + 1) 5 X 10"^ 




K{g ^ 83^1) B{T{1S) ^ fi+fi-){0^^'^\^Si)) = 5 x lO'^ 




K{g ^ 8^Si) B{T{2S) ^ ^i+|I-){oJ^^^\^Sl)) = 1 x lO-^ 




K{g 8^ Si) (Of = (2J + 1) 3 X IQ-^ 


iO 



FIGURES 

FIG. 1. The differential cross section for pp — ipipX versus pt('0) at the Tevatron. 
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